Collective modes of quantum dot ensembles in microcavities 
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Emission spectra of quantum dot arrays in zero-dimensional microcavities are studied theoreti- 
cally, and it is shown that they are determined by the competition between the formation of the 
collective superradiant mode and inhomogeneous broadening. The random sources method for the 
calculation of photoluminescence spectra under a non-resonant pumping is developed, and a mi- 
croscopic justification of the random sources method within a framework of the standard diagram 
technique is given. The emission spectra of a microcavity are analyzed with allowance for the spread 
of exciton states energies caused by an inhomogeneous distribution of quantum dots and a tunneling 
between them. It is demonstrated that in the case of a strong tunneling coupling the luminescence 
spectra are sensitive to the geometric positions of the dots, and the collective mode can, under 
certain conditions, be stabilized by the random tunnel junctions. 

PACS numbers: 42.50.Ct, 42.50.Pq, 78.66.-m, 78.67.Hc 



I. INTRODUCTION 

[ Cavity quantum electrodynamics is one of the most ac- 
.tual directions in the field of modern optics. It studies sys- 
tems where a photon localized in one or several directions 
interacts with elementary excitations of a media P, Q ■ 
'Beams of atoms and single atoms and molecules can stand 
,out for the media. Quantum properties of radiation and 
■matter are most brightly manifest themselfs in these stud- 
ies (see, e.g., Q and references therein). 

Semiconductor nanostructures serve as the solid state 
.analogs of such atomic systems. Investigation of the opti- 
■cal properties of semiconductor nanostructures, and, in 
particular, structures with quantum dots, so-called ar- 
'tificial atoms, is one of the rapidly developing areas of 
|the modern solid state physics 0, H, [1] . The effective in- 
.teraction between localized photon modes and excitons 
"(electron-hole excitations) is attained in semiconductor 
.quantum microcavities . Light-matter interaction being 
multiply enhanced by a cavity allows to efficiently visual- 
|ize eigenstates of nanostructures Q- 

' Lately a remarkable progress in obtaining the strong 
.coupling regime between an exciton localized in a single 
'quantum dot (QD) and a photon trapped in three dimen- 
Isions by a cavity is reached 0, [l^, [HI, • In the case 
■where the damping of a photon mode caused by the imper- 
fections of a microcavity and a non-radiative damping of 
■an exciton become smaller as compared to their coupling 
constant, new eigenstates of the system, zero-dimensional 
exciton polaritons, are formed. Eigenstates structure and 
emission spectra in such cavities have been studied the- 
oretically in a number of works, see e.g. [l^, [ij, [IB]- It 
was assumed in these works that only one quantum dot 
is in the strong coupling regime, the possible presence of 
other emitters has been taken into account phenomeno- 
logically [Tsj . 

In this regard, the study of optical properties of mi- 
crocavities where several emitters with energies near the 
photon resonance are placed is of apparent interest. The 



present paper is aimed at the theoretical investigation of 
this problem. It can be expected, that the interaction 
with li ght will reinforce definite (symmetric or superradi- 
ant [i3) modes in much the same way as it takes 
place in planar microcavities and in optical cavities with 
cold atoms 0, However, as it will be shown here, the 
interaction of light with localized excitons in semiconduc- 
tor quantum dots has a number of fundamental features 
being most brightly manifested in the presence of the tun- 
neling coupling between the dots. 

The paper structure is as follows. The model of zero- 
dimensional microcavity with a quantum dot array is for- 
mulated in Sec. |TT] and a formalism to calculate the pho- 
toluminescence (PL) spectra of such structures based on 
the random sources method is set forth. This method 
is microscopically justified within the framework of the 
standard diagram technique for non-equilibrium systems. 
Third section of the paper is devoted to the investiga- 
tion of different models of exciton generation in resonant 
and non-resonant quantum dots and related features in 
the luminescence spectra. The calculated emission spectra 
of the disordered quantum dot arrays with allowance for 
the dispersion of localized excitons resonance energies are 
given in the forth section. Photoluminescence of quantum 
dot arrays with tunneling coupling is investigated theo- 
retically in Sec. [V] Main results of the work are briefly 
discussed in Sec. IVll 



II. MODEL 

Consider N quantum dots, placed in a zero-dimensional 
microcavity. It is convenient to write the system Hamilto- 
nian in the secondary quantization representation by in- 
troducing the operators c (c^), describing the annihilation 
(creation) of a localized photon, and operators (a|), 
describing annihilation (creation) of excitons, localized in 
quantum dots (i = 1 . . . JV): 
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FIG. 1: Photoluminescence formation scheme in the micro- 
cavity with quantum dots (QD). Subscript i denotes lowest in 
energy exciton states in dots. A non-resonant pumping (pump) 
to the excited "dark" states, denoted by a multi-index fe, not 
coupled with the photon mode (cavity mode) is assumed. Ra- 
diation escapes microcavity due to the tunneling of a photon 
mode through the mirrors. 



Here Wphot is the cavity photon mode energy, coi are the 
exciton energies in quantum dots, the Hermitian matrix 
tij {tij = t* ^) describes the tunneling coupling between 
the quantum dots, and constants Vi describe an interac- 
tion of localized excitons with light. 

It is assumed that excitons tunnel between the dots as a 
whole because Coulomb interaction prevents charge sepa- 
ration 

BHi^ ■ Such a situation can be realized in large 
quantum dots, whose typical scales are larger as compared 
to the exciton Bohr radius. Independent tunneling of an 
electron and a hole leads to an additional decay channel 
of the exciton resonance and can be taken into account by 
standard methods. 

Let us assume for simplicity that the wavefunctions of 
exciton states can be chosen as real. Thus, by matching 
the signs of the wavefunctions, the quantities Vi can be 
made positive. Note, that such a sign change determines 
the signs of the constants tij. Possible difference of the 
signs of the tunnelling constants and its effect on the opti- 
cal spectra of quantum dots in a microcavity is discussed 
in Sec. El 

The cavity mirrors arc not ideal, therefore photon mode 
turns out to be coupled with the continuum of states prop- 
agating outside of the cavity. On the one hand, it leads 
to a finite lifetime of a photon inside the cavity, and on 
the other, to the possibility to detect the emission of the 
system. 

Schematic illustration of the photoluminescence forma- 
tion in the considered system is given in Fig. [TJ It is 
assumed that a non-resonant excitation takes place, i.e. 
the excitons or free electron-hole pairs are generated in 
the excited states whose energies arc far from the pho- 
ton mode. Pumping leads to the formation of the steady 
distribution function of carriers in the excited states in 
quantum dots. These states are "dark" , because they are 
far in energy from the photon mode. The relaxation of 
carriers and excitons from the excited states to the lowest 
ones being in the resonance with photon mode takes place 
due to the interaction with phonons. The emission caused 
by the photon escape away from the cavity is detected. 



A. Random sources method 

In order to calculate photoluminescence spectrum of 
quantum dot ensemble we have extended the random 
sources method developed in Refs. [U, for the quan- 
tum well structures. It consists of the determination of 
the electric field induced as a result of exciton generation 
in quantum dots, characterized by the amplitudes of the 
random sources Fi . Photoluminescence spectrum is deter- 
mined by an intensity of the field in a cavity mode aver- 
aged over the realizations of random sources. Therefore, 
within the framework of this method the inhomogeneous 
system of equations 



N 



(wphot -uj- i7phot)£^ = ViPj 



(2) 



i=l 



describing oscillations at a frequency lo of coupled oscil- 
lators, photon and excitons, under the action of driving 
forces Fi is solved. Here, the quantity E is proportional to 
the field amplitude in the photon mode of the microcavity, 
and the quantities ViPi are proportional to the coordinate 
integrated contributions of excitons in iih quantum dot 
to the media polarization. Eqs. ^ can be obtained from 
the Heisenberg equations of motion for the operators c, 
aj, whose dynamics is governed by the Hamiltonian ([T]), 
and non-radiative dampings of exciton 7^, photon damp- 
ing 7phot and amplitudes Fi are added into Eq. ([2]) in a 
phenomenological way. 

It is convenient to introduce the eigenfrequencies f2„i, 
determined from the homogeneous system ([2]), as well 

as the eigenvectors [cj"^\, cf'' . . . C^"^] = [E,Pi... Pat], 
which represent Hopfield coefficients Q for exciton polari- 
tons in a microcavity. The electric field in a photon mode 
is then given by 

N+l r^{m) 



E{oj) = 



Here the scalar product (C(™'|F) _ 
defined. Photoluminescence intensity is given by the ex 
pression 



I{u)^T\E{u)\^ = T Y 



where the correlation matrix S equals to 



(3) 



(4) 



and double angular brackets denote averaging over the 
random sources realizations. Factor T is proportional to 
the mirrors transmission coefficient and links the intensity 
of the field inside and outside of the cavity. Expression ^ 
shows that the luminescence intensity has, as a function 
of frequency u, poles at the eigenfrequencies 17*^ and il™/, 

and is proportional to the photonic fraction C^^^C^J^, 
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as well as to the random sources correlation matrix S. The 
calculation of quantities Sij is only possible by using the 
microscopic approach, which is developed in the following 
Sec. Ill BI It is shown below that in the relevant case of ex- 
citon quantization as a whole and for the not very closely 
positioned quantum dots matrix S is diagonal, Sij oc Sij. 

We note that the formulae ^ for the photolumines- 
cence spectrum can be simplified in the strong coupling 
regime where the imaginary parts of eigenfrequencies 
Tm = Imfim are negligible as compared with the dis- 
tances between the real parts, i.e. 



(to ^ to'). (5) 



Thus, in the sum in Eq. ([3]) one can neglect the terms with 
to 7^ to', which leads to the expression 



N+l 



(™) 1 2 



(CM|5|CM) 



m— 1 ^ ' '"' 



_ (6) 
generalizing the results of Ref. [2^ to the case of a quan- 
tum dot system placed in a zero-dimensional microcav- 
ity. Condition (O means that the polaritons are well 
defined and, therefore, the structure of terms in Eq. ^ 
can be interpreted in a following way: first two factors 
describe photon emission spectrum from wth polariton 
state (product of a photonic fraction in this state by a 
Lorentzian describing the density of states of a given po- 
lariton mode), and a third factor describes a steady pop- 
ulation of this state, being proportional to the ratio of 
the polariton generation rate (C'™-' jS'lC'™-') and polari- 
ton damping Tm- 



B. Green's function method 

Microscopic description of photoluminescence of quan- 
tum dot arrays in microcavities is possible by using the 
L.V. Keldysh diagram technique. In order to shorten the 
notation we neglect the tunneling coupling between quan- 
tum dots, and in the end of the present section we com- 
ment on exciton tunneling influence. 

One can readily check that the emission intensity is con- 
nected with the Green's function of the photon outside of 
the cavity in a following way 



-hup lim ( — 



(7) 



where 7 is a fictitious damping needed to introduce the 
steady distribution function of photons outside the cavity. 

Under low intensity non-resonant pumping the micro- 
cavity emission is proportional to the population of the 
excited exciton states. Therefore, in Dyson equation for 
Green's functions D~~, one can totally neglect the 



pumping, consequently. 



1 



17 



D++ = -(D — )* 



(8) 



and in equation for -Dg J it is enough to take into account 
only linear in excited exciton states populations contribu- 
tions. As a result, the sought for Green's function can be 




FIG. 2: Diagram representation of Eq. ((9}. Dashed lines de- 
note plioton Green's function outside of the cavity, wavy lines 
its Green's function inside the cavity, solid lines are exciton 
Green's functions in quantum dots calculated without interac- 
tion with light, dotted lines are phonon's Green's functions. 



recast as a sum of two diagrams in Fig. [5] (see also [22 
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NQ + G-^+^^^il + Na) . (9) 



Here the quantities gq describe the link between the cavity 
mode and a continuum of photon states outside of the cav- 
ity, V is Green's function of the photon mode inside the 
cavity, calculated with the allowance for the light-matter 
interaction, G~~ {Gf^) are Green's functions of excitons 
in ith quantum dot, is Green's function of excitons 

in the excited states which are numbered by a subscript fe 
(in a general case it runs through the values belonging to 
both discrete and continuous spectrum). Matrix elements 
U^f^ {U^_^^) in Eq. © describe the transitions from the 
state i to the state k (and vice versa), accompanied by 
an emission or absorption of a phonon with the wavevec- 
tor Q and frequency fig, Nq is the phonon distribution 
function. Phonon damping in Eq. ^ is neglected. 

Exciton Green's functions entering Eq. ([9]) and calcu- 
lated in the lowest in the pumping power approximation 
have form 



Gi ,,, — 



1 



U! - LUi+ 17i 



-7—1 Gf^ — -(Gj^, 



(10) 



Gfc+ = -27r/fc<5(^-£fc), 

where fk is the distribution function in the excited states 
formed as a result of the non-resonant pumping, and Sk 
is the spectrum of the excited states. A simple shape of 
Green's functions is related with the fact that the exciton 
population in the state i at weak pumping is negligible, 
the damping of the excited exciton states is neglected also. 

One has to determine photon Green's function in- 
side the microcavity self-consistently, with allowance for 
both light-matter interaction, described by the Hamilto- 
nian ([1]) , and for the damping of the photon mode caused 
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by photon tunnelling through the mirrors. As it has been 
done above, wc disregard the populations of the excited 
exciton states and of the photon mode. By solving the 
corresponding Dyson equations we have 

vz- 



(11) 



Here 7phot(w) = J2q \9q\'^H^ ~ '-^q) the photon mode 
damping caused by the escape of photon outside of the 
microcavity, small renormalization of the mode frequency 
due to the coupling with an environment is neglected. 

Substituting Eq. © into Eq. ([7]), integrating over fre- 
quency and passing to the limit one can arrive to 
Eq. ([3]). At the same time factor T in Eq. ([3]) turns out 
to be proportional to hujq\gq\'^. Elements of the random 
sources correlation matrix S [Eq. up to a common 
factor equal to 

k,Q 

{l + NQ)5{u + ojQ-ojk)]. (12) 

Therefore, according to its microscopic meaning, Sij {uj) is 
the exciton generation matrix (i.e. the rate of its density 
matrix change) in quantum dots. It is worth noting that 
here we do not issue the challenge of the precise quantita- 
tive description of the carrier relaxation in the quantum 
dot arrays. It is enough to assume that expression ([T^ 
describes the last phonon-assistcd transition in a quan- 
tum dot before the formation of an exciton in the ground 
state [Hi. 

Inasmuch as the expression (C'-'"' |S'|C'-'" -'), determin- 
ing in accordance with Eq. ([3]) photoluminescence spec- 
trum, constitutes a quadratic form with respect to the 
Hopfield coefficients c\"^\ and Eq. is bilinear with 
respect to U^^, ^k-^i^ formulae for the generation 
matrix Sij keeps its shape in an arbitrary basis of exci- 
ton states i, j. Exciton states energies do not enter into 
the expression for Sij at all, therefore, it can be shown 
that Eq. (fT^ holds with the allowance for the tunneling 
coupling between quantum dots as well. 

Let us compare diagonal and off-diagonal elements of 
the matrix 5*. In the relevant case of the exciton localiza- 
tion as a whole the typical distances between the centers 
of quantum dots d exceed exciton Bohr radius. One can 
check that in these conditions Sij <^ Sa {i 7^ j) provided 
the wavelength of the phonon which causes the transitions 
from the excited state k to the states i and j is smaller 
than d. Estimations show that for the acoustic phonon en- 
ergy hujQ ~ 10 meV (which corresponds to the typical en- 
ergy gaps in quantum dots) phonon wavelength amounts 
to the order of 2 nm, i.e. it is significantly smaller than 
the typical distances between centres of the neighboring 
quantum dots. It means that the off-diagonal elements of 
matrix Sij Eq. are small as compared to the diagonal 
ones, thus 



S,^ 



S i Sij , 



(13) 



and the exciton generation in different quantum dots can 
be considered as independent. 
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FIG. 3: Photoluminescence spectra of the quantum dot system 
in a microcavity. Solid curve (1) was calculated for A*" = 1, 
ti^phot = uji,V = 80 /leV, 71 = 17 fieV, 7phot = 50 ^eV, dashed 
curve (2) was calculated for N — 2, o^phot =1^1, 1^2 — = 
0.5 meV, 72 — 71, Si — S2, dotted curve (3) corresponds 
to 72 = 71, Si =0,5*2 7^ and dash-dotted curve (4) to 
72 = 150/ieV, and Si = 0, 5*2 / (unspecified parameters 
are the same as for the first curve). The energy of the photon 
mode is used as the reference point, the spectra are normalized 
so that the maximum values coincide. The tunneling between 



the dots is neglected, ti^- 



0. 



III. NONRESONANT QUANTUM DOTS 
INFLUENCE ON THE PHOTOLUMINESCENCE 
SPECTRA 

The photoluminescence spectrum ^ is determined by 
the structure of cxciton-polariton states and the rate of 
polariton generation in these states. In the present sec- 
tion we analyze how the emission spectra of the cavity 
with a single quantum dot in resonance with photon mode 
is affected by the presence of other, non-resonant, quan- 
tum dots. Here and in the following sections we take the 
exciton-photon coupling constants identical for all quan- 
tum dots, Vi = V > 0, and assume the generation matrix 
Sij to be diagonal. 

The quantum dot photoluminescence in the microcav- 
ity was theoretically studied in a recent work p^ . In the 
case where only one dot interacts effectively with the cav- 
ity mode, the results of Ref. [l^ lead in the limit of weak 
excitation to Eq. ([3]). An influence of nonrcsonant quan- 
tum dots, with \uji — cijphotl 3> V, is taken into account in 
Ref. [l5| only phenomenologically: generation and follow- 
ing emission of excitons in these dots are reduced to the 
effective pumping to the photon mode, which, in turn, in- 
teracts only with the resonant quantum dots. The pump- 
ing to the photon mode can be formally described by 
adding a random source term i^phot to the right-hand side 
of the first of Eqs. the photoluminescence spectrum 
for = 1 is then given by expression 



I{uj) oc 



S, 



phot I 



V^Si 



(14) 



where Sphot ((^phof^p*hot))- advantage of our ap- 
proach is the possibility to take into account explicitly 



5 



all quantum dots interacting with the cavity mode. As 
a result, the quantity S'phot can be related to the exciton 
generation rate, and the contribution of other dots, lying 
beyond the scope of Eq. can be found. 

Fig.[3]shows the photoluminescence spectra correspond- 
ing to the different mechanisms of exciton generation. The 
calculation was carried out for the typical values of the 
quantum microcavity parameters |12l |. indicated in the 
caption to the Figure. Curve 1 was calculated for the 
case where only a single quantum dot is in the micro- 
cavity, so that the spectra demonstrates a doublet corre- 
sponding to the strong-coupling regime. The energies of 
non-interacting photon and exciton modes are assumed to 
be identical and chosen as a reference point, the pumping 
takes place only to the exciton state. 

The curve 2 demonstrates the photoluminescence spec- 
trum for the system where a second dot, detuned at the 
value UJ2 — o^phot = 0.5 meV ^ V = 80 ^eV from the 
cavity resonance, was added, the exciton generation rate 
is the same in both dots, 5*2 = Si. Tunneling interac- 
tion between the dots is neglected. One can see from the 
comparison of curves 1 and 2, that in the frequency range 
\uj — WphotI ^ |w2 — WphotI the presence of a second quan- 
tum dot leads only to a small distortion of the spectrum 
explained by the alteration of polariton eigenfrequencies. 
Analysis shows that the relative magnitudes of the peaks 
in this doublet strongly depend on the detuning uji — Wphot 
and on the ratio of the generation rates Si/ 82- The rise 
of the extra peak in the spectra at the exciton resonant 
frequency in the second dot UJ2 is more important. 

This peak at the frequency 102 becomes dominant in the 
spectrum (curve 3) when the generation of excitons takes 
place only in the second dot, i.e. Si = 0, S2 0. Since 
the second dot acts also as a photon emitter, the doublet 
near the frequency tui remains in the spectrum, however, 
its intensity is smaller as compared to that of the peak at 
the frequency uj2- The shape of the doublet is described 
by Eq. HH) with Si = and 



S'phot — 5*2 



\U! -UJ2+ 172 P ' 



(15) 



so that S'phot <^ 'S'2 when |a;— a;2+i72| ^ V. Note that the 
depth of the dip between the peaks in the doublet for the 
curve 3 is considerably larger than for the first two curves. 
This distinctive feature of pumping "to the photon mode" 
is related to the factor \ui — uji + i'jxl'^ in the numerator 
of Eq. (|14p , which is minimal at the frequency uj = lui. 

At last, the curve 4 is obtained for = 0, S'2 7^ and 
72 ^71- This calculation models the situation when 
the generation of excitons takes place, for example, in a 
continuum of the states (wetting layer) and a direct gen- 
eration of excitons in the resonant state of first quantum 
dot is inefficient. Comparing curves 3 and 4 one can see, 
that the competition between the emission from resonant 
and nonresonant exciton states exists at Si = 0. 

Strong-coupling regime holds for the parameters used: 
two polariton frequencies are Re f2± uji ±V with the 
decay rates r± ~ 71 -I- 7phot < V. The third "polari- 
ton" frequency is obviously close to the frequency of the 
nonresonant quantum dot. The luminescence spectrum is 
well described by Eq. ^ , in this case the structure of this 



expression is so that the peak values of the spectra at the 
frequencies Kc ilm are proportional to l/T^^. The values 
of PL integrated over the frequency round each peak are 
proportional to l/Tm, since the stationary number of the 
polaritons in the state m is the ratio of their generation 
and decay rates. Therefore the emission from the m-th 
polariton state is more efficient for smaller decay rates 
Tm- Consequently, the direct emission from the state 2 is 
more efficient if the quantity 72 is small, and the growth 
of 72 leads to the increasing emission from the resonant 
state 1, in agreement with Fig. [3] 

Thus, we have demonstrated that the luminescence 
spectrum of the quantum dots in a strong-coupling regime 
considerably changes when additional, non-resonant exci- 
ton states exist in the system. The features at the fre- 
quency (jJ2 rise in the spectrum against a background of 
the resonant doublet near the frequency uJi , and the ratio 
of the photoluminescence intensities at the frequencies toi 
and ll)2 strongly depends on the parameters of the prob- 
lem. 



IV. ROLE OF INHOMOGENEOUS 
BROADENING 

Now we proceed to the case where several quantum dots 
with the energies close to the photon mode energy are 
present in the microcavity. In the given section we ana- 
lyze the role of the exciton resonant frequencies dispersion 
due to the inhomogeneous broadening in the ensemble, 
neglecting the tunneling coupling between the dots, the 
influence of the tunneling is studied in the next section. 

It is instructive for further consideration to recall the 
results (isj for identical quantum dots without tunneling, 
i.e. UJi = ojx, li = Ix and tij = 0. One can easily show 
that the luminescence spectrum for such a system is still 
described by Eq. with S2 = 0. There are two mixed 
(polariton) states with the frequencies 



Lox + t^phot - i(7x + 7phot) , 
^± = 7i =t 



i^x - t^phot - i(7x - 7phot) 



(16) 



in this case. All other — 1 exciton states do not in- 
teract with the photon mode and are not manifested in 
the luminescence spectrum. It is evident from Eq. (|16p 
that the splitting between the energies of polariton states 
grows proportionally to ^/N . 

This result can be interpreted by analyzing the eigen- 



vectors P*^*^ = 



Pi 



N 



for the ensemble of not inter- 



acting with light quantum dots, determined by the Hamil- 
tonian ([T|) with V = Q. The energy spectrum is A^-fold 
degenerate for oji = lox and t^.j = 0, so all the eigenfre- 
quencies are equal to t^x, and the eigenvectors P^*^ can 
be chosen as follows: 



P, 



(1) 



1 



N 



J 



.N: 



0, i>2 



(17) 



It follows from the first of Eqs. that for Vi = V the 
strength of the photon interaction with quantum dots is 
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FIG. 4: Photoluminescence spectra for A*' = 9 QDs with inho- 
mogeneous broadening taken into account. Panels (a), (b), (c) 
were calculated for the cavity mode frequency (.jphot = a) (cD is 
the average value of the exciton resonant frequencies) and inho- 
mogeneous broadening a^, = 0, a^^ — 40 ^teV, a^^ — 200 /leV, 
respectively. Vertical lines indicate the positions of exciton 
energies when light-matter interaction is neglected. Other pa- 
rameters are as follows: ^/NV — 80 fieV, 7phot ~ 50 fJ.eV, 
"fi = 17 fj.eV. The pump is assumed to be identical for all 
QDs. Vertical scale is the same for all panels, the spectra are 
normalized so that their maximum values coincide. 



characterized by the quantity X]j=i -^ji i-^- maxi- 
mum when excitons in different dots oscillate in the same 
phase. Therefore, the light interacts only with one of 
the mode s t he symmetrical superradiant mode P^^^ 
H, 0, [M m, ill- This interaction leads to the 

formation of the exciton-polaritons with eigenfrequencies 
(fTB|) , and all other — 1 exciton states are optically in- 
active. 

The light-matter interaction can be conveniently de- 
scribed in the basis of the states (|17p with allowance for 
the dispersion of the quantum dot resonant frequencies. It 
is evident that the polariton eigenstates are determined by 
the competition of the exciton interaction with the cavity 
mode and the inhomogeneous broadening. The first ef- 
fect keeps the states in the form ([TT]) , and the second one 
splits the symmetrical state to the individual localized ex- 
citon modes. The two peak structure of the luminescence 
spectrum should be destroyed when the dispersion of the 
resonant frequencies is comparable to \/NV rather than 
to V. 

This competition is illustrated by the Fig. [H where 
three photoluminescence spectra calculated for N = 9 dif- 
ferent quantum dots with tunneling neglected are shown. 
The panel (a) corresponds to zero inhomogeneous broad- 
ening, so the spectrum consists of two peaks splitted by 



2\/NV. Nine exciton frequencies were randomly cho- 
sen from the Gaussian distribution with the dispersion 
o'uj = 40 /icV for the panel (b). The dispersion of the fre- 
quencies is sufficiently large, it is greater than the light- 
matter coupling constant for the individual quantum dot, 
but smaller, than VNV, = 40 //eV > V « 27 /zeV, 
however, the two peak structure of the spectrum is still 
well resolved, with the depth of the dip between the peaks 
being of the same order as for the panel (a). Calculation 
demonstrates that the anticrossing of the modes happens 
when the detuning of the microcavity is changed. For 
the last panel (c) of the Fig. 2] inhomogeneous broadening 
is large, = 200 /xeV > VNV = 80 fj,eV. The separate 
peaks corresponding to individual dots in a weak-coupling 
regime with the cavity mode are resolved here, the alter- 
ation of the detuning does not lead to an anticrossing. 

Note, that the transition between the weak- and strong- 
coupling regimes in the quantum dot ensembles with large 
homogeneous broadening (decay rate) 7phot,7i > is 
also determined by the ratio of the decay parameters and 
the quantity Consequently, if the characteris- 

tic broadenings of the photon mode and the quantum dot 
spectra are smaller than \/NV, the microcavity is in the 
strong-coupling regime and the emission spectrum has a 
doublet structure. Complex multi-peak spectrum corre- 
sponding to the isolated dots in a weak coupling regime 
takes place in the opposite case. 

V. THE ROLE OF THE EXCITON TUNNELING 
BETWEEN THE QUANTUM DOTS 

In this section the influence of the exciton tunneling be- 
tween the quantum dots on the photoluminescence spec- 
tra is analyzed. In the general case the tunneling coupling 
between the dots can obviously lead to the complete al- 
teration of the eigenmodes Eq. ^T7\ . and to the complex, 
multi-peak photoluminescence spectrum, similar to that 
shown on Fig. ID^c). However, it turns out that the shape 
of the spectrum strongly depends on the specific structure 
of the matrix ti^, i.e. on the character of the tunneling 
couplings between the dots. Two possible limiting cases 
are studied below: random tunneling only between neigh- 
bor dots and random tunneling between all pairs of quan- 
tum dots. Note that the sign of the tunneling contants 
tij is determined by the form of the wavefunctions corre- 
sponding to the exciton states in the dots i and j. If these 
states are identical and of lowest energy then < 0, in 
the general case the sign of can be arbitrary. 

The shape of the microcavity luminescence spectrum 
is qualitatively determined by the distortion degree of the 
superradiant (fully symmetric) mode (|17p , as has been as- 
certained in the previous section. For the following analy- 
sis it is convenient to choose the eigenvectors P^'^ in ([TT]) 
with i > 2, corresponding to optically inactive polariton 
states, as 

^ ^(-^.M - 5,,,), z = 2...N, j = l...N, (18) 

and introduce the tunneling Hamiltonian matrix elements 
{i\t\j) in the basis (HH), {i\t\j) = EL=i ^/''^i,™^™' ■ 
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FIG. 5: Photoluminescence spectra for A'' = 9 QDs, tunneling 
between the dots is taken into account. Curves (a) and (b) are 
calculated for tunneling possible between each pair of QDs, 
t = —80 /ieV t = 0, respectively, and curve () is calculated 
for tunneling only between nearest neighbors, i = —80 fieV. 
Tunneling constants dispersion is at = 80 /leV, other param- 
eters are the same as for Fig. |4]. Vertical lines indicate the 
positions of exciton energies when light-matter interaction is 
neglected. Curve (a) corresponds to tjphot ~ Si (so the cavity 
is tuned to the tunnel-split-off mode frequency); curves (b) and 
(c) correspond to Uphot = i^i- The dispersion of exciton reso- 
nant frequencies uji is neglected. Vertical scale is the same for 
all panels, the spectra are normalized so that their maximum 
values coincide. 



shape, see Fig[5l^a). The curve shown there, is calculated 
for the cavity mode being tuned to the tunneling-split-off 
exciton state. Besides, the spectrum also possesses a fine 
structure in the energy region of the remaining exciton 
states, not shown on the Figure. Note that for the calcu- 
lation presented the tunneling constants dispersion is of 
the same order as the photon coupling constant with the 
superradiant mode, at = VnV. When N increases, the 
strong-coupling regime becomes more stable because both 
the efFective coupling constant and the distance between 
the spit-ofF mode and the other spectrum grow. 

The structure of the spectrum dramatically changes at 
i = (tunneling coupling constants between the dots be- 
ing of the different signs) or when the tunneling takes 
place only between the nearest neighbors (non-zero and 
random elements in the matrix tij are only below and 
above the diagonal). One can readily show that all the 
matrix elements are then of the same order and the 

situation is similar to the inhomogeneous dispersion of the 
energies of the exciton in individual quantum dots. The 
spectra has a doublet structure when at < VnV, in the 
opposite case the eigenmodes are completely altered, in 
agreement with Fig. [5{b) (tunneling constants of differ- 
ent signs, mean value t = 0), (c) (tunneling in the linear 
geometry between nearest neighbors). 

The results presented are relevant for the case, where 
not very large number of quantum dots is present in the 
volume of the microcavity, ^ 1 ... 10. When N is 
large and the dots are chaotically distributed, so that the 
quantities ti j, tOi are random, the structure of the po- 
lariton eigenstates becomes the same as for the planar 
microcavity with inhomogeneous broadening taken into 
account 0, [23|. In case the dots form an ordered two- 
dimensional superlattice, the optical properties of the sys- 
tem are determined by the structure of the Bloch states 
of an exciton 1281 1. 



One can show that if all the random quantities tij = tj^i 
are independent and identically distributed, the averaged 
over disorder realizations values of are equal to 

IJm^iN-l)!, W\3)^0, (19) 

and all the other matrix elements are in order of t with 
root-mean-square deviations being equal to at up to the 
factors of the order of unity. Here i and at are the mean 
value and the dispersion of the random quantity tij, re- 
spectively. Such a situation is realized where the quantum 
dot array is sufhciently dense, so that each dot is coupled 
with all (or at least with most part of the other) dots in 
the ensemble. 

It follows from Eqs. (fT9|) that for sufficiently large N 
and f 7^ 0; the spectrum of the states found with tunnel- 
ing taken into account consists of a single level with the 
energy ei, corresponding to the superradiant mode P(^\ 
split-off by the distance in order of N\t\ from the other 
spectrum. Since the split-off mode weakly interacts with 
others, all N exciton states can still be classified as one, in- 
teracting with the cavity mode and iV — 1 practically opti- 
cally inactive. Therefore, the spectrum keeps the doublet 



VI. CONCLUSIONS 

A theory of the photoluminescence of quantum dot en- 
sembles placed in zero-dimensional microcavities is devel- 
oped in this paper. The random source method, allow- 
ing a transparent description of the photoluminescence of 
such structures as a function of the energy parameters of 
quantum dots and photon mode, is set forth. The ran- 
dom sources approach is microscopically justified within 
the framework of the standard diagram technique for non- 
equilibrium processes, general expressions for the correla- 
tion matrix of the random sources are obtained. 

The influence of the excitation parameters, cavity mode 
energy position and the energy spectrum of the quantum 
dots on the emission spectra is studied. The luminescence 
spectra arc analyzed for the quantum microcavity where 
one dot is in the strong coupling regime with the photon 
mode and energy of all the other ones are far from the cav- 
ity resonance. In this regime the so-called non-resonant 
pumping to the photon mode, introduced in Ref. [lB| for 
the description of the experiments @ is justified both 
qualitatively and quantitatively. 

The structure of the eigenmodes and emission spectra 
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in case where all the dots from ensemble are in resonance 
with the cavity mode is studied in details. An influence 
of the quantum dot exciton energies dispersion and of the 
random tunneling coupling between the dots on the pho- 
toluminescence spectra is analyzed. It is shown, that the 
spectra are determined by the competition of the forma- 
tion of fully symmetric (superradiant) mode and the inho- 
mogeneous broadening in the ensemble. The superradiant 
mode formation is highly sensitive to the geometry of the 
tunncling-coupled quantum dot array, particularly, it is 



strongly suppressed when the tunneling is possible only 
between the nearest neighbors. On the other hand, in 
the dense enough QD ensembles so, that each dot is effec- 
tively coupled with many others, the superradiant mode 
becomes more stable. 
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